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Abstract. RacoonWW is an event generator for e + e~ — > WW — > 4 fermions(+7) 
that includes full tree-level predictions for e + e~ — ► 4/ and e + e~ — > 4/ + 7 as well as 
0(a) corrections to e + e~ — > 4/ in the so-called double-pole approximation. We briefly 
sketch the concept of the calculation on which this generator is based and present some 
numerical results. 



1. Introduction 

At LEP2 and future e + e~ linear colliders, the most important processes to study the 
properties of the W boson are e + e~ — > WW — > Af. For integrated quantities the 
accuracy typically reaches the order of 1% at LEP2 and will even exceed the per-cent 
level at future colliders. To account for this precision in predictions is a non-trivial task. 

High-precision calculations for four-fermion production are complicated for various 
reasons. At the aimed accuracy of some 0.1%, a pure on-shell approximation for 
the W bosons is not acceptable, i.e. the W bosons have to be treated as resonances. 
Since the description of resonances necessarily goes beyond a fixed-order calculation in 
perturbation theory, problems with gauge invariance occur. Discussions of this issue can 
be found in Refs. ||,[§. A second complication arises from the need to take into account 
electroweak radiative corrections beyond the universal corrections. The full treatment 
of the processes e + e~ — > 4/ at the one-loop level is of enormous complexity and involves 
severe theoretical problems with gauge invariance; up to now such results have not been 
published. 

Since lowest-order calculations for e + e~ — ► 4/ have already been extensively 
discussed in the literature (see e.g. Ref. 0] and references therein), we concentrate on 
radiative corrections in the following. Recent results of a full lowest-order calculation 
for real photon emission, e + e~ — > 4/ + 7 are reviewed in Section ^. In Section ^ we 
summarize the strategy and present some results of the event generator RacoonWW, 
which combines full tree-level predictions for e + e~ — > 4/ and e + e~ — > 4/ + 7 with an 
approximation for the virtual 0(a) corrections to four-fermion production via a resonant 
W-boson pair, e + e" — > WW — > 4f. 
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83.4(3) 


6 98(5) 


457(6) 
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running width 


224.6(4) 


84.2(3) 


19.2(1) 


368(6) 


complex mass 


223.9(4) 


83.3(3) 


6.98(5) 


0.460(6) 


e + e~ — > ude"i/ e 7 


constant width 


230.0(4) 


136.5(5) 


84.0(7) 


16.8(5) 


running width 


230.6(4) 


137.3(5) 


95.7(7) 


379(6) 


complex mass 


229.9(4) 


136.4(5) 


84.1(6) 


16.8(5) 



Table 1. Comparison of different finite- width schemes (taken from Ref. 0]) 



2. Full tree-level predictions for e+e — > 4f + 7 

The processes e + e~ — > 4/ + 7 are interesting mainly for two reasons. On the one hand, 
they are a source of (real) 0(a) corrections to gauge-boson pair production with a 
four-fermion final state. On the other hand, they are sensitive to anomalous quartic 
gauge-boson couplings, such as 77WW, 7ZWW, and 77ZZ. In the following we briefly 
summarize some results of Ref. ||, where the building block of the event generator 
Racoon WW |||7j is described that calculates cross sections for e + e~ — > 4/ + 7 with 
arbitrary massless fermions. 

In this event generator different schemes for treating gauge-boson widths are 
implemented. A comparison of results obtained in these different schemes is useful 
in order to get information about the size of gauge-invariance-breaking effects, which 
are present in some finite-width schemes. Table [I] contains some results on the total 
cross section for two semi-leptonic four-fermion final states and a photon, evaluated 
with different finite-width treatments. Similar to the case without photon emission, the 
SU(2)-breaking effects induced by a running width render the predictions totally wrong 
in the TeV range. For a constant width such effects are suppressed, as can be seen from a 
comparison with the results of the complex-mass scheme, which exactly preserves gauge 
invariance. 

Figure [l] shows the photon-energy spectra for some typical four-fermion final states 
that correspond to WW7 production. Apart from the usual soft-photon pole, the spectra 
contain several threshold and peaking structures that are caused by photon emission 
from the initial state. The two relevant classes of diagrams are illustrated in Figure 
Diagrams with the structure of Figure ||a correspond to triple-gauge-boson-production 
subprocesses and yield dominant contributions as long as the two virtual gauge bosons 
can become simultaneously resonant. For instance, WW7 production is dominant for 
£ 7 < 26.3 GeV (224 GeV) for a centre-of-mass (CM) energy of 189 GeV (500 GeV). 
The diagrams of Figure ||b correspond to 7Z production with a subsequent four-particle 
decay of the resonant Z boson mediated by a soft photon or gluon V3. Owing to the 
two-particle kinematics of 7Z production such contributions lead to peaking structures 
around a fixed value of _E 7 , which is located at 72.5 GeV (242 GeV) for a CM energy of 
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Figure 1. Photon-energy spectra (dcr/cLE 7 )/( fb/ GeV) for several processes (taken 
from Ref. g) 




Figure 2. Diagrams for important subprocesses in 4/+7 production (Vi, V% — W, Z, 7, 
^3 = 7jg) 

189 GeV (500 GeV). 

Table [l] and Figure [l| also illustrate the effect of background diagrams, since 
final states that are related by the interchange of muons and electrons differ only by 
background diagrams. While the impact of background diagrams is of the order of 
some per cent for CM energies around 200 GeV, there is a large effect of background 
contributions already at 500 GeV. The main effect is due to forward-scattered e^, which 
is familiar from the results on e + e~ — > Af. More numerical results for e + e~ — > 4/ + 7 
can be found in Refs. 

3. Electroweak radiative corrections to e+e _ — > WW — > 4f 

3.1. Relevance of electroweak corrections 

In the past, Monte Carlo generators for off-shell W-pair production (see e.g. Ref. Q) 
typically included only universal electroweak 0(a) corrections such as the running of the 
electromagnetic coupling, a(q 2 ), leading corrections entering via the p-parameter, the 
Coulomb singularity, which is important near threshold, and mass-singular logarithms 
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a\n{m 2 /Q 2 ) from initial-state radiation, where Q 2 is not determined and has to be 
set to a typical scale for the process under consideration. The size of the remaining 
0(a) contributions is estimated by inspecting on-shell W-pair production, for which 



the exact 0{a) correction and the leading contributions were given in Refs. [[| and flO 
respectively. The difference <5iba — & between an "improved Born approximation" Sjba, 
which is based on the above-mentioned universal corrections, and the corresponding 
full 0(a) correction 5 corresponds to the non-leading corrections and has already been 
discussed in Refs. [@,|TT[. For total cross sections this difference amounts to ~ 1-2% for 
LEP2 energies, but to ~ 10-20% in the TeV range; for distributions the difference is 
even larger in general. Thus, in view of a desired accuracy of some 0.1%, the inclusion 
of non-leading corrections is indispensable. 

3.2. Electroweak corrections in double-pole approximation 

Fortunately, the full off-shell calculation for the processes e + e~ — > WW — > 4/ in 0(a) 
is not needed for most applications. Sufficiently far above the W-pair threshold a good 
approximation can be obtained by taking into account only those contributions that 
are enhanced by two resonant W bosons. The uncertainty from neglecting corrections 
to background diagrams can be estimated to some 0.1%, at least in the absence of 
special enhancement effects such as forward-scattered e^. Doubly- resonant corrections 
to e + e~ — > WW — > 4/ can be classified into two types |2], [12], [L3| : factorizable and 
non-factorizable corrections. 

Factorizable corrections are those that correspond either to W-pair production or to 
W decay. We first focus on virtual factorizable corrections, which are represented by the 
schematic diagram of Figure [3], in which the shaded blobs contain all one-loop corrections 
to the production and decay processes, and the open blobs include the corrections to 
the W propagators. The corresponding matrix element is of the form 

ka _ R+-(k 2 _, k 2 J) R + (k\, k 2 _) k 2 _) 2 2 

{k\ - M*,)(kl - M 2 ,^ kl-M 2 , + kl-M 2 , - - 



doubly-resonant singly-resonant 



non-resonant 
(1) 
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Figure 4. Examples of virtual and real non-factorizable corrections to e + e 
WW -> 4/ 



and the double-pole approximation (DPA) amounts to the replacement 



M 



{k% -M^ + iM w T w )(k 2 _ - M? 

where the originally gauge- dependent numerator R 
independent residue R 



w + iM, 



wFw) 

2 



(2) 



(k%, kt) is replaced by the gauge- 
The one-loop corrections to this residue 
can be deduced from the known results for the pair production H] and the decay |T3| of 
on-shell W bosons. 

The formulation of a consistent DPA for the real corrections, and thus a splitting 
into factorizable and non-factorizable parts, is also possible, but non-trivial. The main 
complication originates from the emission of photons from the resonant W bosons. 
A diagram with a radiating W boson involves two propagators with momenta that 
differ only by the momentum of the emitted photon. If the photon momentum is large 
(E y 3> Tw), the resonances of these two propagators are well separated in phase space, 
and their contributions can be associated with photon radiation from exactly one of 
the production or decay subprocesses. For soft photons (E^ <C I\y), the resonances 
coincide, and the DPA is identical to the one without photon. However, for E^ ~ Tw 
the two resonance factors for the radiating W boson overlap. Although it is possible to 
decompose these factors into contributions associated with the subprocesses, a reliable 
estimate of the accuracy of the corresponding DPA is not obvious. 

Non-factorizable corrections comprise all those doubly-resonant corrections that are 
not yet contained in the factorizable ones, i.e. they include all diagrams involving particle 
exchange between the subprocesses. These corrections do not contain the product of two 
independent Breit-Wigner-type resonances for the W bosons, i.e. the production and 
decay subprocesses are not independent in this case. Simple power-counting arguments 
reveal that such diagrams only lead to doubly-resonant contributions if the exchanged 
particle is a photon with energy E^ < I\v; all other non-factorizable diagrams are 
negligible in DPA. Two relevant diagrams are shown in Figure |], where the full blobs 
represent tree-level subgraphs. We note that diagrams involving photon exchange 
between the W bosons contribute both to factorizable and non-factorizable corrections; 
otherwise the splitting into those parts would not be gauge- invariant. The calculation 
of non-factorizable corrections to e + e~ 



WW — > 4/ was discussed in Refs. [16-18 



in 



detail. A numerical discussion of the sum of virtual and real non-factorizable corrections 
can be found in Refs. 1T7|,|T8 . 
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3.3. Results for O (a) corrections in double-pole approximation 



Different versions of DPA have already been used in the literature |T9|-|22||. For instance, 
following a semi-analytical approach, a consistent application of the DPA for the virtual 
and real corrections to four-lepton production was presented in Ref. |19|]. In such an 
approach, however, it is not possible to apply all experimentally relevant phase-space 
cuts, and effects of recombining photons with nearly collinear charged fermions cannot 
be treated realistically. For the latter reason, in Ref. [HJ the invariant masses of the 



W bosons were defined strictly in terms of the invariant masses of the corresponding 
decay fermion pairs. As a result, the invariant-mass distributions in M± = \fk^ 
received large corrections from final-state radiation, namely —20 MeV, —39 MeV, and 
— 77MeV for r + z/ r , /U + ^ M , and e + z/ e final states at y/s = 184 GeV, respectively. These 
results have been qualitatively confirmed by YFSWW in Ref. pil|| , where the 0{ot) 



corrections to W-pair production |2(J were supplemented by final-state radiation in 



a leading-log approach. Note, however, that the large shifts are due to mass-singular 
logarithms like a ln(mj/Mw), which occur because no recombination of the fermions with 
collinear photons is performed. More realistic definitions of k\, which have to include 
photon recombination, effectively replace the mass-singular logarithms by logarithms 
of a minimum opening angle for collinear photon emission. This expectation was also 



confirmed in Ref. 21 



At present, RacoonWW P,0] is the only Monte Carlo event generator with a 
complete implementation of the 0(a) corrections to e + e~ — > WW — > 4/ in DPA. In 
RacoonWW, only the virtual corrections are treated in DPA, whereas the real-photonic 
corrections are based on the full lowest-order calculation for e + e~ — > 4/ + 7 described 
in Ref. ||. This means, in particular, that only the virtual corrections are split into 
factorizable and non-factorizable corrections. Note that this approach requires a careful 
treatment of IR and mass singularities, since the singularity structure in the virtual 
and real corrections are related to different lowest-order cross sections. For the virtual 
corrections the DPA Born cross section for e + e~ — > WW — > 4/ is relevant, while for the 
real corrections the full Born cross section for e + e~ — > 4/ applies. In RacoonWW one 
can choose between two different methods for the matching between virtual and real 
corrections; one method is based on phase-space slicing, the other on the subtraction 



method described in Ref. [2^]. More details about the RacoonWW approach can be 
found in Refs. @,[F|. 

Figures |5| and ^| show some results of RacoonWW for the semileptonic reaction 
e + e~ — > z/^/x+du at the typical LEP2 CM energy 200 GeV, where 6*w is the W production 
angle (with the W momentum being defined by the total momentum of and and 
Md u is the invariant mass of the quark pair. The precise definition of the input and of 
the photon-recombination procedure, as well as a detailed discussion of more results, 
can be found in Ref. |J. For the photon recombination we first determine the lowest 
invariant mass M 1 f built by an emitted photon and a charged final-state fermion. If 
M 1 f is smaller than a fixed recombination cut M rec the photon momentum is added to 
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Figure 5. Production-angle distribution (dcr/dcos#w)/fb (left) and relative 
correction d/% (right) for e + e~ — > J/^/U+du at -^/s = 200 GeV (taken from Ref. |J) 




Figure 6. Invariant-mass distribution (dcr/dMdu)/(fb/ GeV) of the quark pair (left) 
and relative correction S/% (right) for e + c~ — > i^^t+du at y/s = 200 GeV (taken from 
Ref. §) 



the one of the corresponding fermion /. The curves denoted by "5 GeV" and "25 GeV" 
correspond to the respective values of M rec . While the pro duct ion- angle distribution 
is not very sensitive to M rec , the invariant- mass distribution strongly depends on 
the recombination procedure, as expected from the above discussion. For different 
recombination procedures the maxima of the line shapes differ by up to 30-40 MeV |J. 
As can be seen from Figure || there is a tendency to shift the maxima to larger invariant 
masses if more and more photons are recombined. 
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